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Recent cosmological observations reveal that we are living in a flat accelerated expanding universe. In this 
work we have investigated the nature of the potential compatible with the power law expansion of the universe 
in a self interacting Brans Dicke cosmology with a perfect fluid background and have analyzed whether this 
potential supports the accelerated expansion. It is found that positive power law potential is relevant in this 
scenario and can drive accelerated expansion for negative Brans Dicke coupling parameter u). The evolution of 
the density perturbation is also analyzed in this scenario and is seen that the model allows growing modes for 
negative u). 
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I. INTRODUCTION 

It has been suggested by a number of recent observations 
[|l|] that merely the baryonic matter and dark matter may not be 
able to account for the total density of the universe. This sug- 
gests that either the universe is open (density less than the crit- 
ical density), or, that there is yet another source of energy den- 
sity which makes the universe spatially flat {iltot — !)■ The 
anisotropy of the cosmic microwave background radiation 
(CMBR) as indicated by the recent results from Boomerang 
[^] strongly favour the second possibility of a flat universe 
which makes the issue of missing energy more important. 
Quintessence has been proposed as that missing energy 
density component that along with the matter and baryonic 
density makes the density parameter equal to 1. 

The luminosity-redshift relation observed for the Type-la 
supernovae [||] strongly suggests that, in the present phase, 
the universe is undergoing an accelerated expansion. This is 
supported also by the recent measurements of CMBR and the 
power spectrum of mass perturbations [^. Of course, the sim- 
plest and the most extreme choice in such a case, is the model 
with the vacuum energy density or the cosmological constant 
A 1^, but other options also exist. Models with quintessence 
and cold dark matter (QCDM) is one interesting possibility. 
Basically, quintessence is a dynamical slowly evolving spa- 
tially inhomogeneous component of energy density with neg- 
ative pressure The energy density associated with a scalar 
field Q slowly moving down its potential V can represent a 
simple example of Quintessence [0]. The potential energy of 
the scalar field should dominate over the kinetic energy of the 
field for slow rolling and thereby making the pressure nega- 
tive (pq — ^Q^ — V). For quintessence, the equation of state 
defined by 7q(= ^) lies between and —1. More precisely 

—0.6 > 7q > —1. The deceleration parameter q {= — -^) 
accounting for the acceleration of the universe also lies in the 
same range and —1. Depending on the form of V{Q), jq 
can be constant, slowly varying, rapidly varying or oscillatory 

This quintessence proposal faces two types of problems [Q]. 
One of these problems, (referred as fine tuning problem), is 
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the smallness of the energy density compared to other typical 
particle physics scales. The other problem known as the cos- 
mic coincidence problem, is that although the missing energy 
density and matter density decrease at different rates as the 
universe expands, it appears that the initial condition has to 
be set so precisely that the two densities become comparable 
today. A special form of quintessence field, called the 'tracker 
field', has been proposed to tackle this problem [^. The good 
point about these models is that they have attractor like solu- 
tions that make the present time behaviour nearly independent 
of the initial conditions. 

A scalar field with a potential dominating over its ki- 
netic energy constitutes the simplest example of quintessence. 
Quite a few models have also been suggested with a mini- 
mally coupled scalar field and different types of potentials. 
A purely exponential potential is one of the widely studied 
cases [nffl. In spite of the other advantages the energy den- 
sity is not enough to make up for the missing part. Inverse 
power law is the other potential ( - ) that has been studied 
extensively for quintessence models, particularly for solving 
the cosmic coincidence problem. Though the problems are 
resolved successfully with this potential, the predicted value 
for 7q is not in good agreement with the observed results. 
In search of proper potentials that would eliminate the prob- 
lems, new types of potentials, like Vb[cosh A(/> — 1]^ [ pi] ] and 
Vq smh{aVkoA(l))^ [|^,p^ have been considered, which have 
asymptotic forms like the inverse power law or exponential 
ones. Different physical considerations have lead to the study 
of other types of the potentials also 1 13 1. Recently Saini et al 
[[ij] have reconstructed the potential in context of general rel- 
ativity and minimally coupled quintessence field from the ex- 
pression of the luminosity distance ^^(z) as function of red- 
shift obtained from the observational data. However, none of 
these are entirely free of problems. Hence, there is still a need 
to identify appropriate potentials to explain current observa- 
tions 0. 

As we have mentioned earlier most of the studies to produce 
variable A have been done with a minimally coupled scalar 
field representing the quintessence field. It has been recently 
shown by Pietro and Demaret JI^ ] that for constant scalar field 
equation of state, which is a good approximation for a tracker 
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field solutions, the field equations and the conservation equa- 
tions strongly constrain the scalar field potential, and most of 
the widely used potential for quintessence, such as inverse 
power law one, exponential or the cosine form, are incom- 
patible with these constraints. The minimally coupled self in- 
teracting models will also be ruled out if the observations pre- 
dict that the missing component of the energy density obeys 
an equation of state p — jp with 7 < — 1(/0 > 0) , and these 
sort of equation of state is in reasonable agreement with dif- 
ferent observations [ 16|. These facts motivates toward a more 
general theory like scalar tensor theory. Scaling attractor so- 
lutions are available in the literature with the exponential [ pj| ] 
and power law [ pjp7| ] potentials in non-minimally coupled 
theories. There are different approach also to find tracking 
solutions in general scalar tensor theories with inverse power 
law potentials [18|. Bertolami et. al. [ p^ have found acceler- 
ating solutions with quadratic potentials in Brans Dicke cos- 
mology. In a similar scenario Ritis et al [ pO| ] found a family of 
tracking solutions. We have also found accelerating solutions 
in Brans Dicke Cosmology with a potential which has a time 
dependent mass squared term which has recently become neg- 
ative [ pTj ] . Faraoni [ p2| ] have studied different potentials with 
a non-minimal coupling term ipR^. Like Saini et al [[l4|], 
Boisseau et al p3| ] have reconstructed the potential from the 
luminosity -redshift relation available from the observations in 
context of scalar tensor theory. 

In this work we investigate the nature of the potential rele- 
vant to the power law expansion of the universe in a self in- 
teracting Brans Dicke Cosmology with a perfect fluid distri- 
bution. The Jordan Brans Dicke (JBD) action with a self 
interacting potential and matter is represented by. 



S = 



(1) 



where w is the Brans Dicke parameter and Cm is the La- 
grangian of the matter field. In this theory i plays the role 
of the gravitational constant. This action also matches with 
the low energy string theory action [ |2^ ] for a; = — 1. It 
is checked whether the accelerated cosmic expansion can be 
driven by this potential. The density perturbation is also stud- 
ied to check the consistency of the structure formation sce- 
nario. In the next section we present the exact solution for 
the field equations and investigate the nature of the solution. 
Section 3 discusses the density perturbation part of the solu- 
tion. In the last section we have discussed about the different 
features of the model. 



II. FIELD EQUATIONS AND SOLUTIONS 



The matter content of the universe is composed of perfect 
fluid 



(3) 



where v^v'^ = —1. We assume that the universe is homoge- 
neous, isotropic and spatially flat. Such a universe is described 
by the FRW line-element, 

ds^ = -dt^ + R^{t)[dr'^ + r'^dO'^ + sin^ Odcj)^] (4) 

The Einstein's field equations and the evolution equation for 
the scalar field are given by. 
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The energy conservation equation, that follows from the 
Bianchi identity gives 



R 



(8) 



Among these four equations only three are independent. But 
as there are five unknowns {R,(t), p,p,V) two assumptions 
can be made to match the number of unknowns with the num- 
ber of independent equations. With this freedom, we choose 
the functional form for the time-evolution of the scale factor 
and the scalar field and find the potential that is compatible 
with this choice. We consider solutions for which the scale 
factor and the field (p evolve as power-law functions of time. 



R = R„[ 1] and 



(9) 



where the subscript refers to the values of the parameters at 
the present epoch and to is the present epoch, i.e the age of 
the universe. 

In order to get a solution with accelerated expansion, the 
deceleration parameter, q, has to be negative. This restricts 
the parameter a to be greater than 1. 

The solutions for p and p can be found to be 



P^Pct^-' 
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We start with the Brans Dicke action along with a self in- 
teracting potential and a matter field 



S = J d^x^i^R - ^0"0a - 1/(0) + £,„] 
(We have chosen the unit SttGq = c = 1.) 



(2) 



where. 
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(13) 
(14) 
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The form of the potential allowed by this ansatz is 



V = Vm 



where, 



Vc = 



[&a{uo + 1) - uul3{2 -(3)- 12a^ 



(15) 



(16) 



From @, (|rT|), ([l|) and it is cleai- that the background 
fluid follows an equation of state of the form p = ^^p where 
7/3 is given by 
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Let us now see how the parameters are constrained. First of 
all we are interested in a scenario in which the universe under- 
goes an accelerated expansion. As we have mentioned earlier 
this immediately constrains a to be greater than 1. Further, 
7b, for perfect fluid, is constrained to be between and 1, 
From equation (p^, this implies 2 — 6a</3<2 — 3a. We 
may note that since a > 1, /3 is always negative. 

To clearly specify the nature of the expansion, and the miss- 
ing energy we investigate further the energy density and the 
pressure of the Brans Dicke field. From the field equations, 
the energy density and the pressure of the scalar field turns 
out to be 
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The expression for the energy density and pressure of the 
Brans Dicke are respectively given by. 
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+ (2 - 3a)[/3(l +uj)- (2a + l)]}^t^-' 



(21) 



In order to ensure the positivity condition for both scalar 
field and matter, we find that lo is constrained as. 



l + 2a 



1 < w < 



1 + 2a 



/3 



«(2 - P) 



The allowed region of in lo) parameter space is shown 
in figure 1 for a — 1.1. Now from ( po| ) and (|2l|), it is evident 
that the scalar field describes an equation of state = 
where the 7^ is given by 




FIG. 1. The allowed region in the (/3,tj) parameter space for 
a = 1.1. 
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The expressions for the density parameters for the matter and 
the scalar field, defined respectively by ri,„ = j^r^ and 

= mh [0, are given by 



and 



2a + /3(l + a)-/32(l + cj) 
a(2 - 13) 
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(23) 
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Using these two expressions 7^ can be recast in terms of Vlj, 
and riw, as 
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3a V \l. 



(25) 



As we have mentioned earlier that the result of the Supernova 
Cosmology Project and the High z Survey Project that reveal 
the so-called acceleration of the universe, favours a universe 
with positive cosmological constant as well [||]. Assuming 
flatness in context of general relativity the best fit for these 
data occurred for Vim — 0.28 and JIa = 0.72. Using these 
best fit values along with the bound on (3, it is found that 7^ 
lies within the range 

-0.87 < 70 < -0.47 

In the following two figures we find the allowed region (black) 
in the (a, j3) parameter space which has a 7^ (—0.6 and — 
0.7 with the matter density range within 0.2 < Vim < 0.4. 
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In Brans Dicke theory, the value of the gravitation constant 
is determined by the value of 



(26) 



The scalar field (p approaches to constant value 0o at the 
present time and consequently the inverse of 0o gives the 
Newtonian constant Gn- The rate of change of G at present 



time is given by 



_/3 



t=ta 



Ho where Hq is the Hubble 



parameter at present time (= -§(^0) = ■^)- For any specific 
value of a and f3 admitting accelerated expansion, we need to 
g < 10^^° per year | 




given by HotQ, should be very close to 1. At the same time 
Kaplinghat et al | ^ ] and others [ ^o| ] have also pointed out 
that for power law cosmologies high redshift data indicates a 
to be « 1. In fact in one of our recent investigation we have 
shown that the best fit value for a with SNla data for power 
law cosmology is approximately 1.25 [^T|]. But Hoto is also 
determined by the expression 
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FIG. 4. Hoto vs P for different u) values and a = 1.1 
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FIG. 2. /3 vs a for different flm values and 7^ = —0.6 




FIG. 3. /3 vs a for different SI™, values and 7,^ = —0.7 

In a recent investigation Pont et al estimated the age 
of the universe to to be 14 ± 2 Gyrs, that means a, which is 
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FIG. 5. Hoto vs (3 for different u) values and a = 1.2 

In figure 4 and 5 we plot Hoto with respect to /3 for 
a = 1.1 and 1.2. For this we consider different values 
of uj. For admissible values for the parameters the curves rep- 
resenting the expression Hoto should intersect the horizontal 
line representing a constant value of a at least once within the 
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allowed range of (3. It is interesting to find that for both the 
graphs the curves intersect the straightline representing a at 
least once when oj lies within the range —2.0 < uj < —1.5. 
in a different investigation Banerjee and Pavon j^3|] have also 
arrived at the same range of uj while looking for late time ac- 
celeration in BD theory without a potential. In fact, in other 
investigations [ pTp3| ] in scalar tensor theories studying late 
time acceleration similar conclusions are made. In all these 
analyses the constraint coming from the solar system mea- 
surements (lli > 600) is not satisfied. We would discuss about 
this point in the discussion section. 



III. DENSITY PERTURBATION 

While the accelerated expansion is supported by the Brans 
Dicke theory with an appropriate form of the potential , it is 
worth checking whether the issue of structure formation is 
modified by the dynamics of the Brans Dicke field treated 
here. We are particularly interested in the evolution of the 
density perturbation in the context of Brans Dicke theory with 
the self interacting potential. Our objective here is to find the 
perturbations of the field equations (5)-(8) and analyze the be- 
haviour of the relevant variables for the accelerated expansion 
of the universe. In reference [jsj] a similar treatment is done 
to analyze the behaviour of the perturbed solutions in the orig- 
inal Brans Dicke theory. 

As our interest is in cosmological density perturbations, we 
will be using the temporal component of the field equations. 
Let us denote by Sg^,y, SRoo and 6T'^'^ perturbations in g^^, 
Rqo and T™, respectively. Denoting Sg^^, by /i^^, SRqq can 



be expressed as 



6R, 
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(28) 



Similarly the perturbations in the trace of the energy- 
momentum tensor can be expressed as. 



T ^T + ST ^T + Sp-35p 



(29) 



The perturbation of the d' Alembertian of the Brans Dicke field 
is given by 
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SD<p = 5<p + RRh''^4> - ^hkk^ 



(30) 



The following parameters are used for the relevant perturba- 
tions 



hkk = R'h 

5(t> = X(f> where A << 1 
6p — Ap where A << 1 



(31) 



where /i(t) ,X{t) and A (i) are the perturbed gravitational 
field, scalar field and matter energy density respectively. 

As was mentioned in the earlier section, only three of the 
four field equations are independent. Hence, three equations 
are used for perturbation. They are 



_ P 37B(l-f a;)-|-2-fa; 

R ^ 6 2io + 3 



1 



2 - m 



2uj + 3 
p(l-37s) 2-m 



where 



2a; + 3 2w + 3 



/j + 3(l+7B)|p = 



2-/3 . 
7s = — 1 



and TO = 



3a 

/3 



(32) 
(33) 
(34) 

(35) 
(36) 



Using the above parameters, the set of perturbed equations 
take the form 



h R, 
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(37) 
V^A 

(38) 
(39) 



where is the comoving velocity of the fluid. The growth of 
perturbations primarily occur in the post-recombination era. 
The pressure acting on the matter is negligible during this era 
and hence, it is appropriate to consider the matter to be pres- 
sureless i.e, 73 = 0. From (43) and (44), a and (3 can be 
expressed in terms of to 



/3- 



and a 



TO 
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3(1 -to) 



(40) 
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Hence, the perturbed equations for the pressureless case are 



h R; p , , 2 + a; 
- + — /i = -(A - A) 
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In order to solve the above equation we assume the following 
form 



(43) 



(44) 



A-A = /(t); f{t)^Ct': A = xt' 

where ^ and x constants. 

With the above substitution we find that 9 — 6 and 



The perturbation in the four-velocity Sv'' is to be set null, 
which can be done by one infinitesimal gauge transformation. 
Let us suppose that the perturbation in the scalar field behave 
as plane waves 



X{x,t) — A(t) exp(— ifc.a;) 



(45) 



where k is the wave number of perturbation. Using equation 
(52) and substituting the solutions (9), (12) and (19) in equa- 
tion (50) and (51) we get 
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Combining these two equations and neglecting the higher 
order terms beyond (as the modes are analyzed in the 
asymptotic region) we get, 
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(50) 



xo^ + o[xiC3 - 1) + aci - C3)] + Ciix -0 + 02^ = 

(51) 

The solution for 9 from the above equation 

± Vixil - Cs) + acs - Ci)}2 - 4x{C4x -0 + C2O] 

(52) 

Bertolami et al [|l9|| have considered the growing modes for 
density perturbation in the asymptotic limit of \lu\ >> 1. The 
solution of 9 in this limit for our case is 



1 - m V X 




1 — m 
2lo 



TO ^ 



1-TOX, 



} 



(53) 



This asymptotic value of 9 for large \uj\ depends only on 
the power of the potential m. Since f3 is always negative, 
m(= ^-^) is a positive number. (Note that m — ^^j^ 7^ 1 as 
both a and (3 blow up for m = 1.) So, 9^ represents the grow- 
ing mode for £, < X- Clearly for m — 2 the mode matches 
with the result given by Bertolami et al [|l9|]. However, we 
found that u; should lie in the range —1.5 and —2.0. So in this 
connection it is worthwhile to check if growing modes exist 
with this limitation. For uj = —1.8 we can calculate 9± from 
the equation (|5^. We find that with this value of uj also 9^ 
represents a growing mode. We, of course, need to make suit- 
able choice for the values of ^, x and m to achieve this. Thus 
in our case perturbations can grow with time for range of jcjl 
allowed in this model. 



IV. DISCUSSION AND CONCLUSIONS 

In this work we have investigated the nature of potential 
relevant to the power law expansion of the universe in a self- 
interacting Brans Dicke Cosmology with a perfect fluid dis- 
tribution. The age of the universe and the time variability of 
the gravitational coupling are also calculated. The value of the 
parameters are constrained from different physical conditions. 
We have graphically represented the permissible values of the 
parameters. We have also studied the evolution of density per- 
turbations in this model. 

From our analysis we draw the following conclusions: 
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1 . Accelerated growth of the scale factor can be driven by 
this positive power law potential with 7^ lying in the 
range —0.47 > 70 > —0.87. This range for 7^ is very 
well consistent with the observations. 

2. The gravitational coupling grows with time which 
agrees quite well with the observational facts [^vj]. 

3. This model also allows growing modes for the energy 
density perturbation of matter implying that the dynam- 
ics of the self interacting Brans Dicke field does not up- 
set the structure formation scenario. 

This model depends entirely on three parameters a, (3 and cj. 
The parameter a is only constrained to be greater than 1 by 
the fact that the universe is accelerating. Depending on the 
value of a, the value of /3 gets restricted within a certain range. 
The value of lo depends on these two parameters. In order to 
ensure the weak (positivity) energy conditions as well as to 
produce the observed values of Vim, a and 7^, we require ui 
to be in the range —2 < u; < —1.5. Although this value of u; 
does not satisfy the solar system bound of {u; > 600) this is a 
generic difficulty in most of the investigations done in the con- 
text of scalar tensor theories [ |3l| , p3| ]. In a completely different 
work in BD theory without potential Banerjee and Pavon 
also had a similar range for uj. The work of Bertolami and 
Martins | |l9| ] is an exception where they have obtained accel- 
erated solution with large |a;|. They have, however, not con- 
sidered the positive energy condition for the energy density of 
the scalar field. A large value of necessitates, violation of 
this condition in these models. The violation of the positive 
energy condition is particularly serious because the contribu- 
tion of the energy density from matter is subdominant with 
respect to the missing energy density as is suggested by ob- 
servation. On the other hand this is not the only case when 
the solar system limit is violated. There are other evidences in 
literature where small value of \uj\ has been supported. The ex- 
tended inflationary model suggested by La and Steinhardt 
requires ui to be 20. In order to explain the structure formation 
scenario successfully in scalar tensor theory, the constraint on 
Lu is not also compatible with the solar system test [|3^]. 

Aother imporant point to note is that as the universe ex- 
pands in a constant power law fashion in this model, it has 
an acceleration also in the radiation dominated era which up- 
sets the primeval nucleosynthesis scenario. This problem also 
occurs if the universe is considered to have a power law ex- 
pansion and no compromise can be done between the bigbang 
nucleosynthesis and present acceleratedphase in general rela- 
tivity as is shown by Kaplinghat et al [E^ . Infact the existing 
supernova observational data provides sufficient evidence for 
the fact that the universe was decelerating in the near past [ |3^ ] . 
In a recent work of Banerjee and Pavon [ p3| ] it has been shown 
that in BD theory the nucleosynthesis problem can be avoided 
by considering w to be a function of cf). But in that case also ui 
asymptotically acquires a small negative value to have a late 
time accelerating phase. 



In the end we can say that though the model presented here 
describes the present day universe quite successfully, it dis- 
agrees at two other important point of observation. While 
the age of the universe, density parameters and 7^ calculated 
here satisfies the observed constraints, the power law acceler- 
ated expansion of the universe upsets the standard decelrating 
phase of the universe until recent past and the important so- 
lar system limit on lu. So to get rid of these two problem 
one approach is to consider oj a function of cf) so that local 
inhomogeneities give rise to large uj and one could get both 
accelerating and decelerating phase. We wish to address this 
problem in our next work. 
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